On completely positive non-Markovian evolution 
of a d-level system 



Andrzej Kossakowski 

Institute of Physics 
Nicholaus Copernicus University 
87-100 Torun, Poland 
kossak@fizyka.umk.pl 

and 

Rolando Rebolledo * 

Laboratorio de Andlisis Estocdtico 
Facultad de Matemdticas 
Pontificia Universidad Catolica de Chile 
Casilla 306, Santiago 22, Chile 
rreboUe&puc.cl 



Abstract 

A sufficient condition for non-Markovian master equation which en- 
sures the complete positivity of the resulting time evolution is presented. 



1 Introduction 

An open system is one coupled to an external environment [UH]. The interaction 
between the system and its environment leads to phenomena of decoherence 
and dissipation, and for this reason recently receive intense consideration in 
quantum information, where decoherence is viewed as a fundamental obstacle 
to the construction of quantum information processors [3| . In principle, the von 
Neumann equation for the total density matrix of the system and the reservoir 
provides complete predictions for all the observables. However, this equation 
is in practice impossible to solve since all degrees of freedom of the reservoir 
have to be taken into account. Main efforts have focused in deducing the time 
evolution of the reduced state density matrix. This is the aim of the well-known 
exact theory of subsystem dynamics due to Nakajima-Zwanzig ([H [S]) which 
relies in a generalized (non-Markovian) master equation approach. 



*Partially supported by PBCT-ACT13 and Dircccion de Relaciones Internacionales-PUC 



1 



The Nakajima-Zwanzig projection operator method makes possible to derive 
an exact equation for the reduced density from the von Neumann equation of 
the composed system. The resulting generalized master equation -an integrod- 
ifferential equation- is mostly of formal interest since such an exact equation 
can almost never be even written down explicitely in the closed form. In con- 
trast, when one makes the Markovian approximation, i.e., when one neglects the 
reservoir memory effects, the resulting Markovian master equation [6l [7] takes 
a simple form and the required [5] complete positivity of the resulting time evo- 
lution is maintained. The main goal of the theory of open quantum systems is 
a non-Markovian description of the dynamics which at the same time include 
reservoir memory effects and retain complete positivity. 

A variety of non-Markovian master equations have been proposed (cf. [H 
l[T0l[TTl[Tl[T3[Tl[T3[M[T3[Tl[Ta[20l[2B 

[30]). However, the complete positivity of the resulting time evolution is still an 
important problem to be investigated. 

In the present paper a sufficient condition for non-Markovian master equa- 
tions is given, which ensures that the resulting time evolution is completely 
positive. It is shown that this condition is rather difficult to verify in practice. 
The main reason for that is related to the normalization condition of the time 
evolution. This difficulty can be overcomed if one looks first for completely 
positive unnormalized solutions to non-Markovian master equations, while the 
normalization is imposed separately. 

2 Notations 

Let be a d-dimensional Hilbert space with the scalar product (•, •) and ele- 
ments e, X, y,z, . . .. 

The C*-algebra of linear operators on will be denoted by Ma- Elements 
of Md will be denoted by a,b,c, . . . and the unit of Md is Id- The Md is the 
Hilbert space under the scalar product (a, b) = tr {a*b). 

The C*-algebra of linear maps from Md into Md will be denoted by £{Md), 
its elements are A, B,C, . . . and the identity map in £,{Md) will be denoted by 
id. The conjugation (duality) •'^ in £,{Md) is defined by the relation: 

{A*a,b) ^ {a,Ab), (2.1) 

for all a, 6 G Md- 

This operation endows the following property: the relations 

Aid = Id, Lld^O, (2.2) 

and 

tr (A#a) = tr (a) , tr (L#a) = 0, (2.3) 

are equivalent. 

The cone of all completely positive maps on Md will be denoted by B'^{Md). 
Finally, if At € £,{Md), t > 0, then the Laplace transform of At will be 
denoted by Ap. 
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3 Non-Markovian master equations 



The reduced dynamics can be studied equivalently in the Schrodinger or the 
Heinsenberg pictures. Suppose that At : Md Md describes the reduced 
dynamics in the Heisenberg picture, then it should satisfy the foUowing condi- 
tions: At € B^{Md), Atld = Id, for all t > 0, and Aq = limj^o At = id. In the 
Schrodinger picture these relations are given in terms of Af , i > 0. 

In the present section, the reduced dynamics is investigated under the as- 
sumption that At is the solution of a non-Markovian master equation of the 
form: 

rf^t . /■* . , . 

LAt+ dsLt-sAs, (3.1) 



dt 

with the initial condition Aq ~ id, where 

La = i[h,a] + Fa-^{F{ld),a}, (3.2) 

and h = h* G Md, F €E B~^{Md), that is, L is the generator of a completely 
positive semigroup. 

The normalization condition Af ld — Id. implies the equality 

Ltld^O. (3.3) 



A non-Markovian master equation of the form (|3.ip can be easily derived from 
the Heisenberg equation for the composed system by the Nakajima-Zwanzig 
method under the assumption of factorization of the initial state of the composed 
system and the invariance of the initial reservoir state under the reservoir free 
evolution, c.f. [2]. In this case, La = i[h, a] only, with h — h* £ Md- 

Taking the Laplace transform of (|3.ip one finds: 

(id p - L - Lp)Ap = id. (3.4) 
The equality before implies that both relations below: 

Ap = {pid- L- Lpy^ (3.5) 

and 

Ap{id - L - LpY^ ^ id, (3.6) 

hold. 

It follows from (|3.6p that equation (|3.ip can also be written in the form 



^ = AtL+ f' dsAsLt-s, (3.7) 



and consequently, the dual dynamics becomes: 



^^-LAf+ fdsLt^*A*. (3.8) 
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This means that in the case of non-Markovian master equations there is an 
analogy to the Markovian case. 

To find conditions on L and Lt that ensure that the time evolution At 
resulting from (|3.ip is completely positive for all i > is the fundamental 
problem of non-Markovian master equations. The main result of the current 
paper can be summarized in the following theorem. 

Theorem 1 Let us suppose that At is the solution of the equation p.ip . where 
Lt has the form 

Lt = Bt + Zt, (3.9) 
where Bt E B+{Md) for all t > 0, 

Zta = ^^{Bt{ld),a} + i[ht,a], (3.10) 

and ht — h*^, then At is completely positive for all t > if the solution of the 
normalization equation 



dNt , , , 

LNt+ dsZt^sN,, (3.11) 



dt 

with the initial condition Nq — id, is completely positive for all t > 0. 

Proof. It follows from (|3.ip that the Laplace transform Ap of At is given by 
the formula 

Ap^iidp- L- Zp- Bp)-\ (3.12) 
and satisfies the equation 

Ap = [idp- L- ZpY^ + [id p - L - Zp^^ BpAp. (3.13) 
It follows from p.l3p and (|3.1ip that (|3.ip can be written in the form: 

pt /'t — U 

At^Nt+ du dsNt-u-sBuAs. (3.14) 
Jo Jo 

If Nt is completely positive for all t > 0, then iterating (|3.14p it is easy to see 
that At is completely positive as well for all t > 0, since Bt £ B^{Md), provided 
the iteration procedure converges, q 

In order to analyze the problems related to the solution of the normalization 
equation let us consider the non-Markovian master equation of the form 



dA 
dt 



= [ dsk{t ~ s){Bt-s - id)As, 

JQ 



(3.15) 



where Bt e B+[Md) and Bt{ld) = Id for aU t > 0, and k{t) > 0. 
The normalization equation takes the form 



^ = - fdsk{t-s)N,, 
at Jo 



(3.16) 
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with the initial condition Nq = id. 
The solution of (|3.16p has the form 



Nt = fit)id, 



where f{t) satisfies the equation 

dfjt) ^ 
dt 



dsk{t — s)f{s), 



(3.17) 



(3.18) 



and /(O) = 1. 

As a particular case, let us choose k{t) in the Lidar-Shabani form, cf. [28], 

ie.. 



k{t) = K 



(3.19) 



In this case one easily finds 



fit) 



-^1-7- 



\1 + Kt), if 7 = 1, 



cosh KtW 7^ — 1) -I — 7^= sinh{Kt\/ 7^ — 1) 



, if < 7 < 1 
if 7 > 1. 



(3.20) 

It follows from (|3.17p that Nt is completely positive if and only if f{t) > for 
alH > for all t > 0, and (|3^ shows that f{t) > for alH > if and only if 
7 > 1. 

The above example clearly indicates that the structure of non-Markovian 
master equations is much more complicated than the Markovian ones. 



4 Modified non-Markovian master equations 

The time evolution (in the Heisenberg picture) is given by the family of maps 
At : Md ^ Md, t > 0, such that At € B+{Md), for all t > 0, (complete pos- 
itivity condition), At{ld) — Id for all t > 0, (normalization condition) and 
Aq := lim(|o At = id. In section 3 it has been shown that if At satisfies equation 
(|3.ip . then the normalization condition can be imposed with no trouble. Indeed, 
if p.2p . (|3.9p and (|3.10p are satisfied, then the normalization condition is triv- 
ially fulfilled. On the other hand, complete positivity of At leads to complete 
positivity of solutions to the normalization equation (|3.10p which is a very diffi- 
cult problem. However one can circumvent the above difficulty in the following 
manner. Let Vt, t > be the family of complete positive maps on Md such that 
Mmt^o Vt = id. If Vt{ld) > for all t > 0, then the maps At, t > 0, defined as 

Atia) = 14(l,)-i/Vt(a)T4(l,)-i/2, (4.I) 
are completely positive and normalized. 
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Let Vt, t > he the solution of the foUowing modified non-Markovian master 
equation: 

^^PVt+ [ dsBt-sVs, (4.2) 
at Jo 

with the initial condition limj^o = id, where P is a completely positive map 
and Bt G B+{Md) for aU t > 0. 
The resolvent of (|42|), 

Vp ^ (id p - P - Bp)-\ (4.3) 

satisfies the equation 

% = (idp-P)-^ + {idp- P)-^ BpVp, (4.4) 

which is the integral form of (|4.2p . 

Iteration of (|4.4[) yields that Vt is completely positive since exp(iP) and Bt 
are completely positive. 

If the solution of (|4.2p satisfies the condition Vt{ld) > for all t > 0, then 
At, t > 0, defined through (|4.ip gives the correct time evolution, i.e., it is 
completely positive and normalized. 

The above approach contains as a special case the semigroup form of the 
dynamics. Let us consider the equation 

^ = LVt + X^ f^dse^'-^^^V,, (4.5) 

where L is the generator of a completely positive semigroup. One easily finds 
the solution of (|4.5p which is of the form: 

Vt = cosh(Ai)e*^, (4.6) 

and Vt(ld) = IdCosh(At). The corresponding normalized evolution At has the 
form 

At = e*^, (4.7) 

that is, it is a semigroup. 
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